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ABSTRACT 


During  the  Voyager  1  spacecraft  flyby  of  Titan  on  November  12, 
1980,  an  intense  broadband  low  frequency  electric  field  noise  was 

1  i 

observed  in  the  inbound  wake.  This  so-called  sheath  noise'  is  believed 

l\  t 

to  be  generated  by  an  ion  beam-plasma  interaction  between  the  ions  in 
the  corotational  Saturnian  magnetospheric  plasma  and  nearly  stationary 
plasma  created  by  the  ionization  of  Titan's  atmosphere.  The  analysis 
is  based  on  the  information  from  Voyager  1  and  reasonable  assumptions. 
The  results  agree  quite  well  with  the  observation.  The  analysis  shows 
that  the  instability  only  occurs  where  the  density  of  the  corotational 


ions  is  comparable  to  the  density  of  ions  originating  from  Titan's 
atmosphere.  The  growth  rate  is  high  enough  to  generate  the  observed 
noise,  and  the  calculated  and  observed  frequency  ranges  are  in  good 
agreement . 
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I.  INTRODUCTION 


During  the  Voyager  1  spacecraft  flyby  of  Saturn's  moon.  Titan,  on 
November  12,  1980,  the  plasma  wave  Instrument  detected  various  kinds  of 
plasma  wave  emissions.  There  were  upper  hybrid  resonance  emissions, 
radio  emissions  from  Saturn  and  a  broad  band  of  low  frequency  electric 
field  noise  [Gurnett  et  al. ,  1981;  1982],  The  purpose  of  this  paper  Is 
to  investigate  the  origin  of  the  low  frequency  electric  field  noise. 

Titan  is  the  only  moon  in  the  solar  system  known  to  have  a  sub¬ 
stantial  atmosphere.  At  the  time  of  the  Voyager  1  flyby,  Titan  was 
located  within  the  outer  magnetosphere  of  Saturn.  It  has  been  found 
that  Titan  has  no  appreciable  intrinsic  magnetic  field.  Therefore,  its 
atmosphere  interacts  directly  with  the  magnetosphere  of  Saturn  [Ness  et 
al.,  1981].  Hartle  et  al.  [1982]  described  a  model  of  this  inter¬ 
action,  which  is  shown  in  Figure  1.  The  Voyager  1  trajectory,  project¬ 
ed  into  the  Saturnian  equatorial  plane,  is  indicated  by  the  line 
labeled  VI  trajectory  near  the  bottom  of  the  diagram.  The  electric 
field  noise  occurs  in  the  cross-hatched  region  from  about  0532:30  to 
0538:30  Spacecraft-Event-Time  (SCET).  This  noise  was  called  sheath 
noise  by  Gurnett  et  al.,  because  it  is  qualitatively  similar  to  the 
electrostatic  wave  seen  in  the  raagnetosheath  at  Earth  and  in  the  iono- 
sheath  at  Venus  [Gurnett  et  al.,  1982;  Rodriguez,  1979;  Scarf  et  al., 
1980].  A  similar  kind  of  noise  was  also  observed  near  the  upstream 


edge  of  the  ionized  gas  clouds  produced  by  the  AMPTE  (Active  Magneto- 
spheric  Plasma  Tracer  Explorers)  solar  wind  ion  releases  [Gurnett  et 
al.  ,  1985;  1986].  In  these  cases,  the  origin  of  the  noise  was  explain¬ 
ed  by  an  ion  beam-plasma  interaction  between  the  nearly  stationary  ions 
produced  by  the  gas  cloud  and  the  rapidly  flowing  solar  wind  protons. 

In  the  present  case,  the  sheath  noise  is  believed  to  be  generated 
by  an  ion  beam-plasma  interaction  between  the  corotational  Saturnian 
magnetospheric  plasma  and  the  nearly  stationary  plasma  created  by  ioni¬ 
zation  of  Titan's  atmosphere.  In  Section  II  of  this  paper  we  will 
summarize  the  observations  related  to  the  wave  phenomena  and  a  detailed 
analysis  will  be  made  in  Section  III. 


II.  SUMMARY  OF  OBSERVATIONS 


The  properties  of  magnetospheric  plasma  just  outside  of  the  wake 
region  are  fundamental  to  the  magnetospheric  interaction  with  Titan's 
atmosphere.  It  has  been  known  from  the  Voyager  1  experiments  that 
there  are  at  least  two  major  components,  H*  and  N+,  in  the  magneto¬ 
sphere  near  the  orbit  of  Titan.  In  addition,  N2+  and  ^CN3-  were  also 
possibly  present  [Hartle  et  al.,  1982].  Neubauer  et  al.  [1984]  sum¬ 
marized  the  basic  properties  of  the  magnetospheric  plasma  near  the 
orbit  of  Titan.  The  number  density  of  protons  was  about  0.1  cm-3,  and 
of  nitrogen  ions,  N+,  was  about  0.2  cm-3.  The  other  possible  compo¬ 
nents  were  minor  and  will  be  ignored  in  the  analysis  of  this  paper. 

The  temperatures  of  H+  and  N3"  ions  were  210  eV  and  2.9  keV,  respective¬ 
ly.  The  electron  number  density  and  temperature  were  about  0.3  cm-3 
and  200  eV.  To  a  good  approximation  the  plasma  flow  velocity  was  in 
the  corotational  direction  (20°  deflected  toward  Saturn).  The  speed 
ranged  from  about  80  to  150  km/s ,  somewhat  less  than  the  corotational 
speed  near  Titan's  orbit,  which  is  200  km/s. 

Titan's  exosphere  has  been  studied  in  detail  by  Hartle  et  al. 
[1982].  The  major  composition  in  the  exosphere  is  molecular  nitrogen, 
N2 ,  and  atomic  hydrogen,  H.  The  densities  of  N2  and  H  were  observed  to 
be  108  cm-3  and  4  x  104  cm-3  at  the  exobase,  which  is  4000  km  from  the 
center  of  Titan  (see  Figure  1)  [Broadfoot  et  al.,  1981].  Assuming  a 
temperature  of  160°K,  Hartle  et  al.  [1982]  derived  the  radial  density 
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distributions  shown  in  Figure  2.  The  density  scale  height  of  N2  was 
about  150  km,  and  that  of  H  was  about  1500  km. 

The  upper  hybrid  resonance  emissions  observed  in  the  vicinity  of 
Titan  provided  a  plasma  density  profile  along  the  trajectory  of  the 
spacecraft.  The  top  panel  of  Figure  3  shows  the  plasma  density  pro¬ 
file.  This  figure  was  reproduced  from  Gurnett  et  al.  [1982].  The 
inbound  sheath,  the  tail  region  and  outbound  sheath,  if  there  is  one, 
were  marked.  The  density  level  of  magnetospheric  plasma,  M3-  and  H3-, 
which  was  about  0.3  cm-3,  is  indicated  by  the  dashed  line.  As  can  be 
seen  in  the  figure,  the  density  was  enhanced  up  to  40  cm-3  in  the 
vicinity  of  Titan.  Since  the  plasma  in  the  wake  originated  mostly  from 
the  photoionization  of  exospheric  particles,  the  plasma  density  profile 
was  basically  governed  by  the  exosphere  density  distribution.  The 
scale  height  of  the  plasma  density  is  about  500  km,  which  indicates 
that  the  densities  of  H  and  N2  must  be  comparable  in  this  region. 
However,  there  is  no  information  about  the  fractions  of  H+  and  N3-,  or 
any  other  minor  species,  such  as  N2"3"  or  I^CN3". 

The  low  frequency  electric  field  intensities  detected  by  the 
plasma  wave  spectrum  analyzer  are  shown  in  the  bottom  panel  of  Figure 
3.  The  sheath  noise  from  0532:30  to  0538:30  SCET,  is  clearly  distin¬ 
guishable  from  the  noise  in  the  tail  region,  from  0539:30  to  0545:30 
SCET.  The  sheath  noise  has  a  broad  frequency  range,  from  10  to  1000 
Hz.  A  representative  spectrum  of  this  noise,  at  the  point  marked  by  an 
arrow  in  Figure  3,  is  shown  in  Figure  4.  The  spectrum  lias  a  broad  peak 
centered  at  a  frequency  of  about  100  to  200  Hz.  The  broadband  electric 


field  strength  is  about  0.5  mVm-i.  As  can  be  seen  from  Figure  3,  the 
frequency  of  the  sheath  noise  has  a  tendency  to  increase  as  the  space¬ 
craft  approached  Titan.  During  this  time  the  plasma  density  increased 
monotonically  along  the  trajectory.  However,  a  clear  cutoff  is  seen  at 
about  0539:00  SCET,  when  the  density  is  about  15  cm'3  to  20  cm-3.  The 
intense  noise  was  only  observed  on  the  inbound  leg,  which  is  on  the 
dayside  of  Titan,  and  was  essentially  absent  on  the  outbound  leg,  which 
is  on  the  nightside. 


III.  ORIGIN  OF  THE  SHEATH  NOISE 


The  plasma  processes  in  the  wake  region  are  known  to  be  complicated 
and  involve  ionization,  mass  loading,  and  ion  acceleration.  To  try  to 
understand  the  origin  of  the  plasma  wave  noise,  a  simple  model  will  be 
used.  Since  the  sheath  noise  occurred  at  the  region  where  the  corota- 
tional  magnetospheric  plasma  flow  streams  through  the  nearly  station¬ 
ary  plasma  from  Titan's  exosphere,  an  ion-beam  plasma  instability  is  a 
likely  candidate  to  generate  the  noise.  A  reduced  (one-dimensional) 
velocity  distribution  for  the  model  is  qualitatively  shown  in  Figure  5. 

Six  groups  of  particles  are  considered:  cold  and  hot  electrons,  magneto- 
spheric  M1-  and  if*",  and  exospheric  N+  and  If*-  from  Titan.  Minor 
constituents  are  neglected.  In  the  assumed  model,  each  component  is 
represented  by  a  drifting  Maxwellian  velocity  distribution.  The  drift 
velocity  of  the  magnetospheric  N*  and  ff1-  is  the  magnetospheric  flow  velo¬ 
city,  Vm,  which  is  taken  to  be  150  km/sec.  The  drift  velocity  of  the 

exospheric  M1-  and  H+  is  approximately  zero.  Because  the  electron  cyclo- 

-►  ■+ 

tron  radius  is  small,  both  the  cold  and  hot  electrons  drift  at  the  E  *  B 
drift  velocity,  Vg  .  Since  there  is  no  information  about  the  drift  velo¬ 
city  of  electrons,  we  assume  that  the  electrons  drift  in  such  a  way  that 
the  total  current  Is  zero.  The  validity  of  this  assumption  can  be 
verified  by  estimating  the  possible  current  density  from  the  variation  of 
the  magnetic  field  over  the  interaction  region.  The  magnetic  field  varies 
about  4  nT  over  4  Titan  radii.  Therefore,  the  current  density  is  not  more 


d 


than  3  *  10-^  A/m“.  Taking  the  particle  number  density  to  be  .3  cm-'*,  it 

can  be  shown  that  the  shift  induced  in  the  relative  velocity  between  the 
electrons  and  ions  is  only  7  km/sec,  which  is  less  than  10%.  The  current 
essentially  does  not  affect  our  instability  analysis.  Before  one  can 
proceed  with  an  instability  analysis,  assumptions  must  be  made  concerning 
several  parameters  that  are  either  not  well  known  or  were  not  measured. 

First,  the  electron  temperature  in  the  interaction  region  is  not  com¬ 
pletely  known.  Most  likely  the  electron  temperature  is  similar  to  that  in 
the  magnetosphere,  except  there  must  be  a  cold  component  associated  with 
the  photoionization  of  exospheric  neutral  particles.  Therefore,  we  assume 
that  there  are  two  groups  of  electrons,  hot  and  cold-  The  hot  electrons 
are  assumed  to  have  a  temperature,  Teh,  of  200  eV,  which  is  the  temperature 
of  the  magnetospher ic  plasma.  The  cold  electrons  are  assumed  to  have  a 
temperature,  Tec,  from  2  eV  to  200  eV.  The  lower  limit  is  the  temperature 
that  would  result  from  photoionization  by  the  ultraviolet  light  and  the 
upper  limit  takes  into  account  the  possible  heating  of  the  photoelectrons 
up  to  magnetospher ic  temperatures.  Finally,  the  temperature  of  the  exo¬ 
spheric  N"*"  and  H"*"  is  unknown.  This  temperature  is  presumably  close  to  the 
exosphere  temperature,  which  is  about  186°K.  The  electron  drift  velocity, 

Vg ,  is  determined  by  the  current  J  =  7.  n.e.V  =  elN  V  -  (N  +  N  )V  ], 

ill  bm  b  sd 

where  and  Ns  are  defined  in  Table  1.  Assuming  J  =  0,  one  obtains 


N,  V 

v  =  _JUL- 

d  N  +N 
b  s 
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A  summary  of  all  parameters  and  their  assumed  values  is  given  in  Table  1. 


A  summary  of  all  parameters  and  their  assumed  values  is  given  in  Table  1. 


Instabilities  associated  with  an  ion  beam-plasma  interaction  were 
first  analyzed  by  Fried  and  Wong  [1966].  An  ion  beam-plasma  interaction 
was  proposed  by  Gurnett  et  al.  [1985]  to  interpret  rhe  low  frequency 
electrostatic  noise  observed  in  an  AMPTE  solar  wind  lithium  release- 
Similarly,  it  can  be  applied  here-  For  a  multispecies  unmagnetized 
plasma  the  electrostatic  instability  is  determined  by  the  reduced  distri¬ 
bution  function, 


rn 

F(v)  =  E  —  F  (v) 
m  s 
b  s 


(2) 


where  ms  and  Fs(v)  are  the  mass  and  distribution  function  of  stfl  species, 
and  the  electron  mass-  The  summation  is  over  all  species.  The  plasma 
is  unstable  when  F(v)  is  sufficiently  double-humped  to  satisfy  the 
penrose  criterion  [Krall  and  Trivelpiece,  1973].  As  can  be  seen  from 
Figure  5,  instabilities  are  in  principle  possible  for  phase  velocities 
between  0  and  Vj  and  between  and  Vra-  However,  the  thermal  velocity  of 
the  beam  ions  is  about  equal  to  Vra,  which  is  larger  than  Vm— .  There¬ 
fore,  the  instability  between  and  Vra  cannot  occur  because  the  beam  ion 
distribution  is  too  broad  to  produce  a  double  hump.  Instability  can  only 
occur  in  association  with  the  relatively  cold  exospheric  ions  at  phase 
velocities  between  0  and  . 

Before  proceeding  with  the  analysis  we  will  consider  the  effect  of 
the  magnetic  field.  In  the  present  problem,  the  magnetic  field  is 
essentially  perpendicular  to  the  plane  of  Titan's  orbit  and  therefore  to 
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the  plasma  flow  direction.  Its  magnitude  Is  about  5  nT.  The  electron 
cyclotron  frequency  is  about  140  Hz.  The  cyclotron  frequency  of  a  proton 
is  less  than  0.1  Hz  and  that  of  an  N*  ion  is  even  less.  The  wave 
frequency  we  consider  here  is  much  larger  than  the  ion  cyclotron 
frequency.  Therefore  the  effect  of  the  magnetic  field  on  the  ion  motions 
is  negligible.  The  electron  cyclotron  frequency  is  almost  in  the  middle 
of  the  wave  frequency  range.  In  general,  electrons  should  be  treated  to 
be  magnetized.  However,  the  gyroradius  of  the  thermal  electron  is  9.5  km 
[Neubauer,  1984].  Even  for  the  photoelectrons,  the  gyroradius  is  an 
order  of  km.  It  is  anticipated  that  the  wavelength  of  the  electrostatic 
noise  will  be  short,  probably  only  several  tens  of  meters.  For  such 
short  wavelengths,  much  less  than  the  cyclotron  radius,  the  wave  motion 
is  basically  not  affected  by  the  magnetic  field.  The  theoretical  proof 
for  this  result  was  given  by  K.  N.  Stepanov  [1959].  Therefore,  we  will 
neglect  the  magnetic  field  in  the  calculation  and  verify  the  consistency 
afterwards.  Another  effect  of  the  magnetic  field  is  that  It  could  affect 
the  distribution  functions  of  the  exospheric  N-1"  and  H+  ions,  which  may 
differ  from  our  assumption.  We  will  discuss  this  effect  later. 

The  linearized  electrostatic  wave  dispersion  relation  for  a 
multispecies  plasma  with  no  magnetic  field  is  given  by 
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where  k  is  the  projection  of  wave  number  on  the  beam  direction,  m  is 


the  complex  frequency  whose  real  part  gives  the  frequency  and  imaginary 


part  gives  the  growth  rate  for  the  instability,  and  ws  is  the  plasma 
frequency  for  sc^  species.  If  Fg(v)  is  a  Maxwellian  distribution,  the 
dispersion  equation  (3)  can  be  written  as 


D(k,u)  =  1  -  j 


V(z  )  - 
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(4) 


where  Xps  is  the  Debye  length  for  sc^  species,  Z'  denotes  the  deriva¬ 
tive  of  the  plasma  dispersion  function  (see  Fried  and  Conte  [1961]  with 
respect  to  its  argument,  and  zs  is 


z 
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(5) 


where  as ,  Vs  are  the  thermal  and  drift  velocity  of  s1-*1  species,  respec¬ 
tively.  The  dispersion  equation  is  solved  numerically  by  using  Salzer's 
method  to  evaluate  the  plasma  dispersion  function  and  Muller's  method 
for  root-finding  [Fried  and  Conte,  1961;  Salzer,  1951;  Muller,  1956], 

Figure  6  shows  a  representative  growth  rate  for  the  ion-beam  plasma 
instability.  The  top  panel  shows  the  wave  number,  k,  normalized  by  the 
Debye  length,  Xp,  as  a  function  of  frequency.  The  bottom  panel  shows 
the  growth  rate,  y,  vs.  frequency.  The  parameters  used  here  are  as 
follows:  Nps  =  Nns,  Ns  =  5  x  Nb  =  1.5  cm-3,  Tec  =  200  eV,  i.e.,  Ngc  = 

0,  Neh  =  Ng  +  Nb>  and  the  values  of  the  parameters  in  Table  I  are  also 
used.  With  the  above  parameters,  the  Debye  length  is  computed  to  be 
about  0.77m.  Taking  kXp  to  be  0.05  or  higher,  the  wavelength  is 


estimated  to  be  about  20m  or  less.  This  wavelength  is  consistent  with 
the  assumption  that  the  magnetic  field  can  be  neglected.  As  seen  in 
Figure  6,  the  growth  rate  is  very  high,  up  to  6  sec-1,  which  means  that 
a  wave  disturbance  can  grow  a  thousand  times  larger  in  1  second.  At 
Voyager  l's  closest  approach,  the  radial  distance  to  the  center  of  Titan 
was  about  6500  km.  As  is  shown  in  Figure  1,  the  unstable  region  in  the 
upstream  side  of  the  spacecraft  path  has  the  scale  size  of  at  least  6500 
km.  The  wave  encountered  by  the  spacecraft  must  start  in  this  region, 
and  propagates  downstream  along  with  the  plasma  flow.  If  a  wave  starts 
1000  km  upstream  of  the  spacecraft  path,  for  example,  it  takes  5  seconds 
(at  a  speed  of  200  km/sec,  including  the  wave  propagation)  to  reach  the 
spacecraft.  During  this  time  the  wave  can  grow  to  1015  times,  which  is 
more  than  enough  to  grow  from  a  thermal  level  to  very  high  amplitudes. 
The  frequency  range,  up  to  200  Hz,  agrees  very  well  with  the  spectrum  in 
Figure  4. 

For  the  calculations  in  Figure  6,  the  wavevector,  k,  is  parallel 
to  the  beam  velocity.  It  can  be  shown  that  for  the  assumed  parameters 
a  wave  propagating  parallel  to  the  beam  velocity  has  the  highest  growth 
rate.  Figure  7  gives  the  maximum  growth  rate,  Ymax>  (maximized  over 
| k | )  as  a  function  of  9,  the  angle  between  k  and  Vm  under  the  condition 
of  Nps  =  Nns>  Ns  =  Njj.  The  solid  curve  is  for  the  case  of  no  cold  elec¬ 
trons,  and  the  dashed  curve  is  for  the  case  Nec  =  Ns.  As  can  be  seen  in 
the  figure,  the  growth  rates  for  both  cases  decrease  monotonically  as  9 
increases,  and  the  instabilities  cease  at  about  85°.  Similar  results 
are  obtained  in  the  variety  of  the  conditions  in  Table  l.  Therefore, 


in  the  rest  of  the  study  we  will,  consider  the  roost  unstable  case,  that 
is ,  9=0. 

In  order  to  study  the  effect  of  the  plasma  density  on  the  insta¬ 
bility,  the  marginal  instability  condition  (y  =  0)  roust  be  investigated 
as  a  function  of  the  plasma-beam  density  ratio,  Ns/Nb.  Figures  8  and  9 
show  the  frequency  of  marginal  instability  as  a  function  of  density 
ratio,  Ng/Njj.  The  dashed  curves  in  these  plots  show  the  frequencies 
corresponding  to  the  maximum  growth  rates.  Figures  10  and  11  show  the 
maximum  growth  rate  versus  Ns/Nb.  The  conditions  for  Figures  8  and  10 
are  that  Nns  =  Nps ,  Tgc  =  Teh  =  200  eV ,  NgC  =  0,  and  Neh  =  Ns  +  Nb. 

The  conditions  for  Figures  9  and  11  are  that  Nns  =  NpS,  Tg*1  =  200  eV, 
Tec  =  2  eV,  Nec  =  Ns,  and  Neh  =  Nb.  The  parameters  not  mentioned  take 
the  values  given  in  Table  1.  These  two  groups  of  figures  cover  the  two 
extreme  limits  for  the  electron  temperatures.  The  two  extremes  give 
very  similar  results.  As  can  be  seen  in  Figures  10  and  11,  the  maximum 
growth  rates  reach  about  10  sec-1,  which  is  quite  high.  The  frequency 
range  extends  up  to  several  hundred  Hz,  which  is  in  good  agreement  witli 
the  observed  frequency  range.  One  of  the  most  important  features  in 
Figures  8  and  9  is  that  the  instability  only  occurs  for  a  limited  range 
of  density  ratio,  Ns/Nb,  between  about  0.001  to  50  for  nec  =  0,  and 
between  about  0.001  to  25  for  Nec  =  Ns.  This  condition  shows  that  the 
instability  cannot  occur  too  far  from  Titan,  where  the  newly  created 
plasma  density  is  too  low  or  too  close  to  Titan  where  the  newly  created 
plasma  density  is  too  high.  This  result  explains  the  region  where  the 
noise  occurs  in  Figure  3.  As  can  be  seen,  the  sheath  noise  started 
about  six  minutes  before  Voyager  I's  closest  approach  to  Titan  where 
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the  density  ratio  is  much  less  than  1,  and  it  ceased  at  the  place  where 
the  density  ratio  is  about  50,  in  good  agreement  with  the  theory.  As 
mentioned  in  Section  II,  the  frequency  range  increased  as  the  density 
increased.  Figures  8  and  9  predict  this  same  trend.  The  unstable  fre¬ 
quencies  start  from  a  few  Hz  when  the  density  ratio  is  about  10-3  and 
increase  monotonically  to  a  few  hundred  Hz  as  the  density  ratio 
increases  to  about  50.  The  effect  of  the  different  fractions  of  N3-  and 
H+  on  the  instability  has  been  studied.  The  results  are  similar  to 
Figures  8  -11  except  that  more  protons  give  a  slightly  larger  frequency 
range. 

No  sheath  noise  was  observed  in  the  outbound  wake.  There  are  two 
reasons  for  this  inbound-outbound  asymmetry.  One  is  a  day-night 
asymmetry.  The  outbound  wake  is  on  the  nightside.  There  is  less 
photoionization  on  the  outbound  wake  than  in  the  inbound  wake,  which  is 
on  the  dayside.  The  other  is  the  direction  of  acceleration  of  the 
newly  created  ions.  This  is  shown  in  Figure  12.  On  the  inbound  side 
of  Titan,  the  ions  are  accelerated  away  from  the  exosphere  and 
therefore  can  escape  downstream.  The  free  escape  caused  the  density 
enhancement.  However,  on  the  outbound  side  of  Titan,  the  ions  are 
accelerated  into  the  exosphere  and  are  absorbed  (see  Gurnett  et  al. 
{1982]  and  Bridge  et  al.  [1981]). 

In  the  above  analysis,  we  have  assumed  that  every  species  of  par¬ 
ticles  has  a  Maxwellian  velocity  distribution.  This  is  certainly  a 
good  assumption  in  the  absence  of  a  static  field.  However,  in  the  pre¬ 
sent  case,  there  is  a  magnetic  field,  which  is  roughly  perpendicular  to 
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the  Titan’s  equatorial  plane,  and  an  electric  field,  which  is  induced 


by  the  corotational  motion  of  the  Saturnian  magnetosphere.  Therefore, 
the  validity  of  this  assumption  should  be  discussed  for  each  species. 

As  is  shown  in  the  upper  panel  of  Figure  13,  if  a  constant  K  field 
is  in  -z  direction  and  an  E  field  is  in  x  direction,  then  the  charged 
particles  will  gyrate  in  the  x-y  plane  in  addition  to  the  drift  motion 
in  y  direction.  In  the  case  of  a  homogeneous  plasma,  the  velocity  dis¬ 


tribution  in  the  x-y  plane  will  be  a  ring  with  the  radius  of  V_ ,  which 
is  equal  to  cE/B,  as  shown  in  the  lower  panel  of  Figure  13.  However, 
if  the  temperature  of  a  particle  is  high,  such  that  the  thermal  velo¬ 
city,  Vth,  is  greater  than  V0,  the  distribution  function  will  be  close 
to  a  single  Maxwellian.  In  the  present  problem,  the  thermal  veloc  ,es 
of  electrons  are  about  850  km/sec  (for  Tec  =  2  eV)  to  8500  km/sec  (for 
Teh  =  200  eV).  The  thermal  velocities  of  the  beam  H+  and  tf4"  are  about 
200  km/sec.  They  are  greater  than  VQ,  which  is  taken  to  be  150  km/sec. 
Therefore,  for  these  four  species,  the  Maxwellian  distribution  is  a 
good  assumption.  Since  the  exospheric  H4"  and  N4"  have  very  low  tempera¬ 
tures,  their  thermal  velocities  are  only  a  few  km/sec,  which  is  much 
less  than  V0.  Therefore,  the  effect  of  the  magnetic  field  on  the 
distribution  of  these  two  species  must  be  considered  in  detail. 

A  ring  distribution  in  the  x-y  plane  can  be  modeled  by  a  function 
as 

f(V  V  V  =  fi(vz)f2(V  V  (6) 
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where  fi(vz)  is  just  a  single  Maxwellian  and 


f2(vx*  V - lA -  exp[-(Vv2  +  (V  -  V  )2-  V  }2/a2]  (7) 

X  y  2  it  V  a  x  y  o  o 

o 


where  a  =  /2kT/m,  and  T,  m  is  the  temperature  and  the  mass  of  the  par¬ 
ticle,  respectively.  A  plot  of  f2(vx,  Vy)  is  shown  in  the  lower  panel 
of  Figure  13.  By  taking  f2(vx,  Vy)  to  be  the  distribution  function  for 
the  exospheric  H+  and  N+,  we  find  that  the  ion  beam-plasma  instability 
is  stabilized. 

To  explain  this,  in  Figure  14  we  plot  the  one-dimensional  distri¬ 
bution  function  F(v)  =  /  f2(Vx,  V)  dVx  for  protons  with  temperature  T  = 
186°K  and  VQ  =  150  km/sec.  The  subscript  y  is  dropped  in  the  one¬ 
dimensional  function.  The  plot  shows  F(v)  for  v  between  -5  km/sec  to 
150  km/sec  because  of  the  symmetry  of  F(v)  about  v  =  V.  This  plot  is 
very  different  from  a  single  Maxwellian  distribution  with  the  same 
temperature  and  VQ  =  0  in  the  following  sense.  Firstly,  the  peak  of 
the  present  distribution  function  near  v  =  0  is  much  lower  than  the 
peak  of  the  single  Maxwellian  because  the  particles  for  the  latter  are 
mostly  concentrated  in  an  interval  of  less  than  10  km/sec,  about  v  =  0, 
while  less  than  10%  of  particles  for  the  former  stay  in  the  same 
interval.  Secondly,  for  the  single  Maxwellian,  the  distribution 
function  is  almost  zero  for  v  >  10  km/sec,  while  for  F(v)  in  Figure  14, 
there  is  tall-like  distribution  with  finite  negative  slope  except  at 
v  =  150  km/sec.  Because  of  this  factor,  no  positive  slope  is  obtained 
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when  we  introduce  an  ion  beam  with  a  temperature  of  210  eV  and  drift 
velocity  of  150  km/sec,  and  electrons  with  temperatures  more  than  2  eV 
and  drift  velocity  0  <  Vq  <  150  km/sec.  Therefore  the  ion  beam-plasma 
instability  cannot  occur  in  this  situation. 

In  doing  this  calculation,  a  homogeneous  plasma  is  assumed.  This 
assumption  is  valid  if  the  scale  length  of  the  interaction  region,  L, 
is  much  larger  than  the  particle  gyroradius,  Rc.  In  the  present  case, 
the  scale  length  is  about  500  km  to  2000  km.  The  lower  limit  is 
obtained  from  the  measured  plasma  density  profile  shown  in  Figure  3, 
and  the  upper  limit  is  basically  the  scale  length  of  hydrogen  radial 
distribution  for  Titan's  exosphere  (see  Figure  2).  The  gyroradius  of 
the  newly  created  exospheric  H+  and  N*  are  248  km  and  3470  km,  respec¬ 
tively.  From  these  estimates  we  see  that  the  ion  gyroradii,  Ke ,  is 
comparable  to  the  scale  length  of  the  interaction  region,  L.  The 
formation  of  a  ring  distribution  requires  that  nearly  the  same  number 
of  particles  move  in  every  direction.  Apparently,  this  requirement  is 
not  met  because  a  large  number  of  particles  would  go  out  of  the  region 
before  completing  one  gyration.  Also,  the  electric  field  acting  on  the 
newly  created  ions  is  not  as  large  as  |Vq  x  Bj  because  of  the  shielding 
effect  of  the  electrons.  From  these  facts  we  believe  that  the  newly 
created  ions  dominate  with  the  velocity  nearly  zero,  and  that  a  single 
Maxwellian  distribution  is  stilL  a  good  assumption. 


IV.  CONCLUSION 


The  low  frequency  electric  field  noise  observed  by  Voyager  1  dur¬ 
ing  the  inbound  crossing  of  Titan's  wake  has  been  analyzed  in  detail. 
The  analysis  shows  that  the  noise  is  generated  by  an  ion  beam-plasma 
instability  caused  by  the  corotational  magnetospheric  ions,  mostly  H+ 
and  NT*",  streaming  through  the  nearly  stationary  plasma  produced  by 
photoionization  of  Titan's  exosphere.  This  mechanism  explains  the 
basic  features  of  the  noise.  By  using  known  (measured)  and  reasonably- 
assumed  parameters,  one  gets  high  growth  rates,  and  good  agreement  with 
the  observed  frequency  range.  The  instability  has  a  limited  range  for 
the  plasma  to  beam  density  ratio,  Ng/N^,,  from  about  0.001  to  50.  This 
explains  why  the  noise  only  occurs  in  the  inbound  wake. 

However,  some  detailed  aspects  of  the  noise  cannot  be  explained 
because  of  the  lack  of  some  information  and  the  complication  of  the 
process  in  the  Titan-magnetospheric  interaction.  For  instance,  the 
calculated  frequency  for  the  maximum  growth  rate  does  not  match  the 
observations  in  detail,  and  also,  the  calculated  frequency  range  for 
the  particular  density  ratio  does  not  follow  observed  density  exactly. 
Actually,  the  detailed  process  in  the  interaction  region  is  not  well 
known.  Also,  for  sufficiently  high  growth  rates,  nonlinear  effects  are 
always  present.  Nonlinear  effects  stabilize  the  system,  saturate  the 
wave  growth  and  broaden  the  frequency  spectrum.  Therefore,  a  simple 


linear  theory  is  not  expected  to  explain  the  sheath  noise  completely. 
Plasma  simulations  are  required  to  fully  investigate  nonlinear  effects 
including  the  role  of  these  waves  in  heating  the  plasma. 

A  similar  kind  of  noise  was  also  observed  in  the  vicinity  of  the 
space  shuttle  [Shawhan  and  Murphy,  1983 J,  apparently  caused  by  an 
interaction  of  gases  from  the  shuttle  with  the  surrounding  ionosphere. 
This  type  of  noise  seems  to  be  a  common  phenomena  when  a  neutral  gas 
interacts  with  a  rapidly  moving  plasma. 
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Quantities 


H+  density  in  beam 
N4-  density  in  beam 
Total  density  in  beam 
Newly  created  H+  density 
Newly  created  N4-  density 
Total  newly  created  density 
H+  temperature  in  beam 
N4"  temperature  in  beam 
Exospheric  H4  temperature 
Exospheric  N4  temperature 
Hot  electron  temperature 
Cold  electron  temperature 
Beam  velocity 
Electron  drift  velocity 
Cold  electron  density 
Hot  electron  density 


FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


The  magnetospheric  plasma  flow  around  Titan  and 
Voyager  1  trajectory,  projected  in  the  Saturnian 
equatorial  plane.  The  electrostatic  noise  was 
observed  in  the  cross-hatched  region.  Since  this 
noise  occurs  in  a  sheath-like  region  around  Titan, 
it  is  called  sheath  noise. 

The  radial  density  distribution  of  Titan's  exosphere, 
normalized  by  the  density  of  H  at  the  exobase,  nQ 
[Hartle  et  al.,  1982]. 

The  electron  density  profile  along  the  Voyager  1 
trajectory  provided  by  UHR  emissions  (top  panel)  and 
the  low  frequency  electric  field  noise  observed  in 
the  vicinity  of  Titan.  The  dashed  line  in  the  top 
panel  indicates  the  total  density  of  the  magneto- 
spheric  plasma.  The  sheath  noise  appeared  from 
0532:30  to  0538:30  SCET,  and  tail  noise  appeared  from 
0539:30  to  0545:30  SCET. 

The  electric  field  spectral  density  of  the  sheath 


noise  at  the  time  marked  by  the  arrow  in  Figure  3. 


Figure  5  A  model  for  analyzing  an  ion  beam-plasma  instability. 

The  magnetospheric  H+  and  fT*"  ions  flow  with  speed  Vm, 
and  electrons  (hot  and  cold)  drift  with  speed  ,  the 
newly  created  ions  are  nearly  stationary.  The  curves 
are  not  drawn  to  scale. 

Figure  6  Plots  of  the  growth  rate,  y,  and  normalized  wave 

number,  k,  as  a  function  of  frequency  under  the  condi¬ 
tions  Ns  =  5  Nfo  =  1.5  cm-^,  Nns  =  NpS,  and  Nec  =  0. 

Figure  7  The  maximum  growth  rates  as  a  function  of  0,  and  the 

angle  between  ic  and  under  the  conditions  Nns  =  Nps 
and  Ns  =  N^.  The  solid  curve  is  for  the  case  Nec  =  0, 
and  the  dashed  curve  is  for  Nec  =  Ns. 

Figure  8  The  frequency  of  marginal  stability,  r  =  0,  as  a  func¬ 

tion  of  the  density  ratio,  Ng/Nj,,  assuming  Nns  =  Nps, 
and  Nec  =  0.  The  dashed  curve  indicates  the  frequency 
associated  with  the  maximum  growth  rate. 

Figure  9  The  same  plot  as  in  Figure  8  with  Nec  =  Ns  and  Tec  = 

2  eV. 


Figure  10 


The  maximum  growth  rate  as  a  function  of  the  density 
ratio,  Ng/Njj,  for  the  case  in  Figure  8. 


Figure  11 


The  maximum  growth  rate  as  a  function  of  the  density 
ratio,  Ns/N{j,  for  the  case  in  Figure  9. 
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